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Abstract. In the context of an inflationary Randall-Sundrum Type II braneworld 
(RS2) we calculate spectral indices and amplitudes of cosmological scalar and tensor 
perturbations, up to second order in slow-roll parameters. Under very simple 
assumptions, extrapolating next-order formulae from first-order calculations in the 
case of a de Sitter brane, we see that the degeneracy between standard and braneworld 
lowest-order consistency equations is broken, thus giving different signatures of early- 
universe inflationary expansion. Using the latest results from WMAP for estimates of 
cosmological observables, it is shown that future data and missions can in principle 
discriminate between standard and braneworld scenarios. 
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1. Introduction 

The idea that the world we live in has more dimensions than we can see dates back to the 
1920s with the works of Kaluza 1J and Klein [2J; however, in order not to violate results 
coming from gravitational experiments, extra dimensions should be compactified and 
very tiny (of order of the Planck scale) and so almost unobservable. During recent years 
new models have been proposed (ADD 0-0, RSI 0, RS2 [7j) which require dimensions 
with compactification scale close to the limit of modern data (< mm, ADD and RSI 
scenarios) or even noncompact dimensions (RS2); in these cases, the visible world is 
confined into a four- dimensional variety (a brane) embedded in a larger spacetime, thus 
called braneworld. Besides regaining classical gravity at low energies, these models have 
many interesting consequences, such as the mass hierarchy problem solution and the 
concrete possibility to test and bound the theory by means of astrophysical (supernovae) 
and cosmological observations (cosmic microwave background temperature fluctuations, 
large-scale structures) and accelerators (high-energy processes available at future LHCs). 
However, the problems opened by this new trend of research are far from being fully 
explored. 

According to the inflationary paradigm, cosmological large-scale structures were 
originated by exponential dilatation of quantum fluctuations up to macroscopic scale, 
thanks to the accelerated expansion driven by a scalar field filling up the early 
universe. The study of microcosm allows us to investigate macrocosm in some sense, so 
cosmological observations are complementary to those with ground-based accelerators. 
Moreover, they clarify the present composition and geometrical structure of the universe, 
as well as those primordial elementary processes constituting the basis of our visible 
world. 

Our aim is to apply the slow-roll formalism to a single field inflationary model on a 
brane in order to compute some observable quantities whose experimental determination 
from CMB data will be greatly improved by the WMAP and Planck missions. In 
particular, we want to compute up to second order scalar and tensor spectral indices 
and perturbation amplitudes in analogy with the four- dimensional standard approach; 
from these observables we will get the fundamental consistency equations which are 
typical of inflationary models. According to the lowest-order results (e.g. [HJ), there 
is no way to discriminate between standard and extra-dimension scenarios, since the 
consistency equation is miraculously the same. The hint of our work is that second- 
order expressions do break this degeneracy, at least formally. The size of the effect 
depends on both physical assumptions and observable quantities, such as the tensor-to- 
scalar-amplitude ratio. Modern observations turn out to be in the range required for 
testing the high-energy braneworld cosmology. 

This paper is organized as follows: in the rest of this section we set up the theoretical 
background; in section 2 we introduce the slow-roll formalism, making some additional 
remarks in section 3; section 4 is devoted to first-order perturbations and in section 5 
we extend the analysis to second order by reasonable guesses; in section 6 observational 
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implications of our results are considered. Section 7 is devoted to discussion. 



1.1. Theoretical background 

We will consider the RS2 scenario [Zj in which matter, and in particular the inflationary 
scalar field, is confined on a single brane while gravity is free to propagate in the whole 
empty bulk. For a complete review of standard inflation, see 9J. The starting point is 
the continuity equation, 

p + 3H(p + p) = 0, (1) 

and the modified Friedmann equations 



H 2 = ^tP(2A + p) + - d , (2) 



6A rv " ' r/ ' a 4 

§ = -£w? + *) + ^ + *)]-!, (3) 

where H is the Hubble parameter restricted to the 3-brane, k 2 = 8n/m1 includes the 
effective Planck mass m 4 ~ 10 19 GeV, A is the brane tension, p is the matter density 
and £ = const, is the dark radiation term which is the time-time component of the 
five-dimensional Weyl tensor projected on the brane. Gravity experiments impose the 
bulk curvature scale to be < 1 mm, that is A 1//4 > 100 GeV. For a homogeneous scalar 
field <j)(t) with potential V we have 

p = ±0 2 + V{4>) , (4) 
p = \tf - V{(j>) . (5) 

In the following we will forget the dark radiation term since during inflation it is strongly 
suppressed by the exponential expansion [T3J [TH] . 

On the brane, the energy-momentum tensor of the scalar field is covariantly 
conserved, = 0, and the field satisfy the same 4D Klein-Gordon equation, 

+ 3H<p + V'((j)) =0. (6) 

A necessary and sufficient condition for inflation to start is a > 0, which is equal to a 
constraint for the equation of state p = wp, 
1 /A + 2p\ 

W < - o > ( 7 ) 



3 V A + p 

which implies 



?<V-^^(^-2V). (8) 



Finally, we note that 



' * 1 2 r 



H = -^( X + P)(P + P) = + 2(A + V)] . (9) 
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2. Slow-roll parameters 

The formalism of slow-roll parameters (hereafter denoted by SR) gives good control 
over the theoretical shape and amplitude of cosmological perturbations. During the 
calculation, formulae containing the smallest power of any of these parameters will be 
referred to as 'first order SR'. In first order SR it is possible to drop the second derivative 
in the equation of motion (JHJ), — —V/3H, the Hubble parameter can be considered 
almost constant and the expansion becomes nearly exponential. Such an approximation 
has been widely exploited in braneworld cosmology [HI EH El and we will call it extreme 
slow roll (ESR). 

In order to refine previous results it is necessary to push onward to second order SR 
and more exactly define the slow- roll parameters. In fact, in the ESR approach one drops 
subdominant dynamical terms thus losing information during the system evolution. 
We define the following slow-roll parameters: 
= H _ 3(X + p)(p + p) 
( H 2 P (p + 2\) ' [W) 

(ID 



V 2 - (12) 




H 2 <p 

Note that e > and, since 

- = # 2 (l-e), (13) 
a 

the condition for inflation is precisely e < 1 and inflation ends exactly when e = 1. 
Moreover, functional dependence of the parameters on and H is the same as in the 
standard case; this is because in both scenarios the definition of inflation is the same, 
a > 0. 

The equation of state for the scalar field becomes 
>(i e -l) + A(|e-l) 



V = P 



P + X 



(14) 



3. Limits 



In the limit in which the brane tension A is much larger than the scalar field energy 
density during inflation, the quadratic term in (J2J) is negligible and one recovers the 
behaviour of the 4D standard cosmological expansion, as one can see by taking the limit 
X/p -> oo in equations Q, ©, 0, ©, © and (JUJ): 

H 2 ^-p, ~^--(p + 3p), 

6 a o 

u><-§, <P 2 <V, (15) 

k 2 - 2 /2 



p_ I _ e -l )p. 
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In the limit p ^> X there is a deviation from the standard equations which, at least in 
principle, can be measured experimentally. 

As has been said before, the slow-roll parameters as functions of H and (f> do not 
change with respect to the standard case; however, since the Friedmann equations are 
different, some relations are no longer valid (for example, the Hamilton- Jacob i equations) 
and one must take care of which expression must be used to define each quantity. 

In the ESR limit, the slow- roll parameters become 

m\ fV'\ 2 AX(V + X) 



16ttV^/ {V + 2X) 2, 
m\V" 2X 

^8^WT2V (17) 
as in [T7]. In the limit 7 C A, equations (JTHJ) and (fTTj) reduce to their standard form 
(see, e.g., PI). 

4. Cosmological perturbations — first order SR 

4 ■ 1 ■ Scalar perturbations 

Taking the definitions of the first-order amplitude and spectral index of the scalar 
perturbations are 

1 H 2 



Svr |0| 



(18) 

k=aH 



1 dln ^I 

n *- ls -dliiP (19) 

where As is evaluated when the wave number k of the perturbation crosses the horizon 
during inflation. Notice that the amplitude has the same H and dependence as in four 
dimensions; this is due to local conservation of the stress-energy tensor on the brane 
[T7] . Defining the number of e-foldings as N((p) = J" ^ H(t')dt', where t e marks the end 
of inflation, k can be written as k(4>) = H(<p)a((f)) = a e H(<p) exp[A r (0)], leading to the 
exact relation 

d In k H , , . 

— (20) 

Thus, to first order the spectral index is the standard one, 

n s -l = 2r]-Ae. (21) 

4-2. Tensor perturbations and the consistency equation 
As regards tensor perturbations, one has ^Sj 

(22) 



k=aH 



2 H „ H 
A T = — —F[- 

din Ai . . 

nrS -55nf- (23) 
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where \x = a/47tA/ (3m|) and 



F(x)= Vl+x 2 - x 2 In ( 1/x + v 7 ! + 1/ 



'X' 



-1/2 



In the limit p <C A, F(x) ~ 1 and one recovers the standard spectrum: 

2 H 



A 



Ut = — 2e . 



The consistency equation is, noting that A\j A\ ~ e, 



ay 



In the limit p ^> A, ~ ^3x/2 and 

6 ^3 \V2 



A, 



\257rm 4 
-3e. 



) (=») 



1/4 



fc=a# 



(24) 
(25) 

(26) 

(27) 
(28) 



As one can see, the tensor index is different from the previous case and the amplitude 
is larger because Kaluza-Klein massive modes contribute in a considerable way in this 
energy regime; however, and rather surprisingly, the consistency equation is the same, 
equation (see |5]). 

The consistency equation reflects the common physical origin of scalar and tensor 
perturbation; it is a typical result from inflation that other theories of structure 
formation are not able to reproduce. 



5. Cosmological perturbations — second order SR 

Assuming that general relativity is valid up to the characteristic scale of the extra 
dimension, several steps have been made towards a full five- dimensional large-scale 
perturbation theory [T4*l fT3| \T§\ I21H |2*T] . An implementation of perturbation theory in 
a brane-inflationary context is still at an early stage. However, it is possible with some 
care to use the well known 4D formalism [22]~|21]- As we want to make a rough estimate 
of braneworld effects on observable anisotropies, we can carry out the calculation with 
little effort by making a somewhat drastic simplification. In particular, we will neglect 
any effects coming from both the total anisotropic stress and the Weyl tensor. The first 
assumption reduces the degrees of freedom of gauge invariant scalar perturbations in 
the longitudinal (conformal Newtonian) gauge The second one closes the system 
of equations allowing us to study brane physics without nonlocal contributions from the 
bulkt. 

f If one chooses a conformally fiat background, ps — Sa~ A = 0, the fluctuation of dark radiation suffers 
an exponential damping during inflation, Sps oc a~ A ~ exp(— 4Ht). 
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This oversimplified background is justified by noting that any improvement of the 
physics would only confirm the breaking of degeneracy between standard and braneworld 
consistency equations, which is our main point. We repeat once again that the 'softness' 
of this breaking and its evidence will strongly depend on the physics, so we cannot say 
hardly anything a priori about its size in more complicated scenarios. However, as a 
last remark, it is important to notice that the effect of this extra-physics (fluctuations in 
Weyl component and anisotropic stress in a high-energy regime) will be more enhanced 
at small scales, since in a conformally and spatially flat background it only mildly 
affects the density perturbations and spectrum at large scales, k <C aH [2TJ |2HJ 
Since inflationary dynamics dominates the large-scale perturbation spectrum, in our 
context we will expect to find results which are close to the 'true' answer coming from 
a complete computation of the full Einstein equations with boundary conditions. 

As we will see soon, the starting point is to find an exact solution of the cosmological 
equations such that the SR parameters are constant, thus ignoring their time derivatives. 
One then perturbs the solution with small variations of the parameters. For example, 
from equation (|rUjl and (fTTj) one gets 

^^-^S' (29) 

r, = H(a, - ( 2 ) , (30) 

infinitesimal if e, rj, £ and ^/H are small enough. Hawkins and Lidsey [27] have found 
four braneworld exact models. Here we will consider just two of them. 
The first one (model 1) has 



a{T) = (T + y/W^iy , (31) 



<j){T) = ~ln(T + VT^T), (32) 

7 

where p > 1/3, j — a/4tt/ (pmf) and T = a/47tA/ (3p 2 ml) (t — t ) is the rescaled time, 
with to being an arbitrary integration constant. Using the notation qr = dq/dT for the 
variable q, we have 

E^= - P . (33) 

The slow-roll parameters are 

1 T 

e = V = — ; , (34) 

e = V^- (35) 



p 2 T 2 - 1 
In the second model (model 2 



a(T) = (AT 2 - If' 2 , (36) 
cosh[ 7 0(T)] = 2T , (37) 
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and 

£ =K i+ iK>- (39) 

14:T 2 ~1 

» = p4tJTI' (40) 
f2 =?( 1 + 4^)(^Tl) ! ' (41) 

In the limit T — > oo, these solutions tend to power-law inflation j2H], with H = p/t 
and e = 77 = £ = l/£>. This occurs both when the brane tension is very large and at 
large times. Actually one can fix the integration constant to such that 0(t = 0) = 0, 
whence to ~ —m^/^/X; with this choice, one can show that the approximation with 
constant SR parameters is valid at late times, for instance when t > 10 3 t4, where 
t4 = mj 1 is the Planck time. Thus we have shown that there exist exact solutions with 
late time constant SR parameters. Anyway, the reader can convince oneself that the 
approximation is good by checking the behaviour of the scale factor a near the origin of 
time (model 1: T ~ 1; model 2: T ~ ^). In both cases, again, one obtains a power law 
(a ~ T p and a ~ t P//2 , respectively) which, combined with the asymptotic behaviour of 
4>(T), generates constant SR parameters. 

5.1. Scalar perturbations 

Let us now see what happens to scalar metric fluctuations. The scalar field is confined 
on the brane and the effective action giving the equation of motion @ reproduces 
standard four-dimensional cosmology. Scalar perturbations, which are related to scalar 
field fluctuations, can thus be treated in the standard way, following step by step the 
procedure we now briefly sketch (e.g. j^]): (1) write the linearly perturbed metric in 
terms of gauge-invariant scalar quantities; (2) compute the effective action of scalar field 
fluctuation and the associated equation of motion; (3) write the perturbation amplitude 
as a function of an exact solution of the equation of motion with constant SR parameters; 
(4) perturb this solution with small variations of the parameters. 

In the extra horizon limit k/aH — > 0, the scalar amplitude for an inflationary model 
with constant SR parameters is j22] 

2-3/2 T(u) / 1 X 1/2-^2 



.4 



s 



5tt T(3/2) \1 



T k=aH 

where v is a combination of SR parameters; it does not matter which of the exact 
models is chosen. Let us now perturb about one of the braneworld solutions with late 
time constant e, assuming that e, rj, £ and (f> / H are small quantities. It can be shown 
that in both the models we are considering one has 4> A /H ~ t -3 at late times. It follows 
that 



A s = [1 - (2C + l)e + Crj\ K 



1 H 



(42) 

k=aH 
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where C ~ —0.73 is a numerical constant. The computation of the spectral index gives 

2A-K d> A 

n s - 1 = [-4e + 277 ~ 8 (C + l)e 2 - 2Cf + 2 (5C + 3)e^] + (2C + 1) — 7^ [1 + O (e, 77)] . (43) 

The first additive term on the right-hand side is the standard expression for the spectral 
index, while the second term comes from the embedding of inflationary cosmology in a 
five-dimensional bulk. For p <C A one has, using (|T5|) . 

~ ^e 2 (44) 
and the extra term can be neglected, thus finding the standard result 

n s - 1 = -4e + 2r] - 8(C + l)e 2 - 2C£ 2 + 2(5C + 3)e^ . (45) 
For p > A, from equations (gj), ©, ©, (EH), (IT4l and (EHjl we get 

12 '* 4 e", (46) 



\m\ H 2 



e = He{e-2r,), (47) 



whence 



n s - 1 = -4e + 2r] - 2(2C + 3)e 2 - 2C^ 2 + 2(5C + 3)er] . (48) 

The only difference with respect to ()45j) is the numerical coefficient in front of e 2 , 
— 8(C + 1) ~ —2.16 in the first case and — 2(2C + 3) ~ —3.08 in the second case. 

5. 2. Tensor perturbations and consistency equations 

Tensor perturbations are more difficult to deal with than scalar ones since gravity is 
free to propagate in the whole five-dimensional spacetime. In general, the graviton zero- 
mode, localized on the brane, interacts with Kaluza-Klein massive modes generating 
an infinite tower of coupled differential equations. In some limits the analysis simplifies 
and the zero-mode decouples; for instance, for a de Sitter brane (p + p = 0) the problem 
can be resolved and one gets equation (j2*2*|) [To] . 

However, some considerations can give semi-quantitative hints about the shape of 
the perturbation amplitude. In fact, we know that at the lowest order in SR parameters 
the amplitude is given by equation f|22[). while in the low-energy limit we must obtain 
the standard 4D amplitude which is known to the next-to-lowest SR order [3*U] : 

At= -4=[l-(C + l)e] — 

5 V^ rn A k=aH 

Our ansatz for tensor amplitude to the second order in SR parameters is then 



2 . . .. . . H ( H 
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The calculation of spectral index and consistency equations is a simple task in the low 
(p <C A) and high (p A) energy limits, labelled by the superscripts (I) and (h), 
respectively. In the first case one gets the standard result [SHI EH] : 

> {0 - - e[2 + 2{2C + 3)e - 4(C + 1)77] (50) 

n s ) . (51) 



IV. 



A 2 
A 2 



A 2 



Since to second order there appear only quantities present in first-order expressions, it 
is reasonable to consider equation (JoTj) not as an extension of ()2b|) but as the consistency 
equation. 

For p ^> A one gets 



n 



(h) 



-e[3 + (2C + 5)e - 4(C + l)rj] . 
Using the approximate relation 



2 j\.rp 

3A 2 



[l-2C(e-T))) 



and the first-order equation ()21|1. we obtain the main result of this work 



n. 



(h) 



2 



3-2-f + (2-C)(l-n s ) 



- 2 



.4 



2 r 

T 



4 



.4 



1 - 0.67^| + 0.91(1 - n s ) 



A 2 



(52) 
(53) 

(54) 
(55) 



The first-order degeneracy between the two regimes is thus (softly) broken. Note that 



Ah) 



to first order n s ' = n s ' = ns- 



We can get a second consistency equation by computing ay = dn^/dlnA; starting 
from equations flSfl), (|52)l. (|25jl. (gTJ) and pp]). For p < A, we have a T ~ 4e(p - e), 



whence, using (j2Tj) and the first-order relation A\jA\ 



ax 



A 2 
' 4 2 



2§ + (»5 " 1)' 



(56) 



When /) > A, (jUJ) gives «t — 3e(2p — e) and so, since A\jA\ ~ 3e/2, one finds (jHBj) 
again. In the light of the results of the previous section, this coincidence is not very 
striking, having used only first-order expressions for the computation. 

A final consistency equation comes from the running of the scalar index, as = 
dns/dlnk. In the low-energy limit one finds as ~ 2(— 4e 2 + 5en — £ 2 ), while in the 
limit p > A we have as — 2(— 2e 2 + 5ep — £ 2 ). These are expressions quadratic in 
the SR parameters, that can be written in terms of observables thanks to the standard 
assumption £ <C max{e, \r)\}, which permits us to neglect the last term, thus obtaining 



a 



a 



(/<) 



2 Arj-i 

3 If 



12^f + 5(n 5 -l) 

32 A.rp 

+ 5(n s -l) 



(57) 
(58) 



For a recent study of the running of the spectral index, see [33] . 
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6. Observational consequences 



Let us now see how much the low and high energy scenarios are different by comparing 
the consistency equations at the two energy scales, equation (joTj) . (jol)j) - (j58j) : 



n, 



(h) 



n 



(0 



a 



(h) 



a 



(0 



a 



(h) 



— a 



(0 



2 A.rp 

34 



\_A\ 
'3A| 



A 5 



+ (l + C)(n 5 -i; 



44 ^ . 



r 



(59) 

(60) 
(61) 



The last equation remains valid if we relax the condition £ <C max{e, by assuming 
instead the parameter £ to be almost constant in the energy scale, ~ fW. Note that 
the running of the tensor index does not give observational evidence of the braneworld. 

The first-year data of the Wilkinson Microwave Anisotropy Probe (WMAP) [3*lj 
[39] confirm the standard scenario of a flat, adiabatic universe with Gaussian, almost 
scale-invariant anisotropies. Bennett et al [S3] put a bound on the amplitudes' ratio, 
r = 16At/A 2 s < 0.90; here we will take the upper limit A\jA\ = 0.06 and the best fit 
for the scalar index of [39:, n$ = 0.93. Substituting these values for the observables in 
equation (joTj) . we have 



a 



a 



(0 



0.01 



(62) 



comparable both with the error in the estimate of Bridle et al jlOj, as = —0.04 ± 0.03, 
coming from the combination of WMAP and 2dFGRS (2 Degree Field Galaxy Redshift 
Survey) data, and with the uncertainty estimate of the Planck mission [H]. However, we 
will get a much smaller effect if the contribution of tensor perturbations is significantly 
suppressed. 

Data analysis carried out for WMAP makes use of the consistency equation 
in order to fix the tensor index and its running in the space of parameters, so a direct 
confrontation between these quantities and an experimental result with an associate 
error is not possible. However, taking the above estimates and putting ny ~ —0.1, 



equation (j59[) gives \(n 



ih) 

T 



71 



(0 



0.02; testing this effect would require an 



experimental uncertainty less than 1% for the tensor index, a very difficult goal to 
hit for the missions of this generation. 

As a last comment we note that the comparison of the observable quantities tit, 
«t and as must be done with the consistency equations, equations ([59* |) - ([5Tjl and not 
through the SR parameters expressions, equations ([HU]). (jB^jl and following. This is 
both because we are dealing with independent expressions and because there is an 
evident ambiguity in relations between SR parameters and observables when considering 
energy-scale finite differences of the quantities of interest. For example, one would have 



„(0 



= 4e 2 , where e would be arbitrarily evaluated along the energy scale. 
Further braneworld effects on the CMB spectrum and confrontation with WMAP 
data are studied e.g. in 
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7. Discussion 

Within the noncompact braneworld Randall-Sundrum model, we have computed the 
main theoretical quantities describing the cosmological perturbations generated from 
a brane-confined inflationary period. In order to find the braneworld signatures in 
cosmological observations we have carried out a semi-quantitative investigation starting 
from previous de Sitter results instead of doing actual 5D calculations, which are very 
complicated due to the presence of the full tower of gravitational modes (that is, the 
brane is not an effectively closed 4D system). In particular, we have written spectral 
indices and consistency equations in the low-energy and high-energy limits using slow- 
roll parameters, to lowest and next-to-lowest order. 

To first order in the SR parameters, the braneworld consistency equation is the same 
as the standard one, equation (J2H|) . and so it is not possible to distinguish between the 
two scenarios. To second order, however, the consistency equations are formally different 
in the two regimes, thus permitting us, in principle, to find observational clues in favour 
of the braneworld model. Even if the present data do not allow such a discrimination, 
new data coming from WMAP and the future Planck satellite could throw more light 
on the problem. 

Finally, there are some generalizations that can give new details of the cosmological 
perturbation structures. In particular, one can of course relax the condition of vanishing 
Weyl or energy-momentum tensor in the bulk |^] , allowing the presence of matter 
outside the brane too. For the case of five-dimensional inflation, see [IH]-(IH]. 
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